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ABSTRACT 



As cells become smaller and, consequently, handoffs more 
frequent, it becomes increasingly important to base the handoff 
decisions on as good information as possible. The purpose or this 
paper is to exploit the information in signal strength 
measurements to improve the quality of the handoff decisions. 

As a vehicle moves, the signal strength changes due to (at 
least) three different phenomena: 

1. Distance from base station. 

2. Slow fading - shadow or lognormal fading. 

3. Fast fading -- Rayleigh fading. 

The fluctuations due to slow and fast fading are usually considered 
as noise in this decision making process. This leads to the 
following tradeoff problem for the averaging interval for the signal 
strength measurements. If the interval is too short, the fading 
fluctuations are not sufficiently smoothed out. If the interval is too 
long, the handoff mechanism is too sluggish. With this tradeoff in 
mind, we present a method to adaptive! y change the averaging 
interval. The method is based on estimating the maximum Doppler 
frequency, the key to the tradeoff. 



/. Introduction 

As cells become smaller and, consequently, handoffs 
more frequent, it becomes increasingly important to base the 
handoff decisions on as good information as possible. The 
purpose of this paper is to exploit the iriformation inherent in 
signal strength measurements to improve the quality of the 
handoff decisions. We seek to obtain useful information from 
signal strength variations which are normally considered only 
as nuisances. 

Asa vehicle moves, the signal strength changes due to (at 
least) three different phenomena: 

1 . Distance from base station. 

2. Slow fading - shadow or lognormal fading. 

3. Fast fading -- Rayleigh fading. 

The decision for handoff is (or should be) primarily 
dependent on the distance dependence (although severe 
shadowing may dictate a handoff). The fluctuations due to 
slow and fast fading are usually considered as noise in this 
decision making process. This leads to the following tradeoff 
problem for the averaging interval for the signal strength 
measurements. If the interval is too short, the fading 



fluctuations are not sufficiently smoothed out If the interval 
is too long, the handoff mechanism is too sluggish. In (1J, 
this tradeoff is quantified We show how to adaprively 
choose the averaging interval. The adaptation uses 
information already available in the signal strength 
measurements in a simple manner. 

The time scale for the Rayleigh fading depends on the 
maximum Doppler frequency (which is proportional to 
vehicle velocity). It would be expected that the time scale for 
the shadow fading should also depend on vehicle velocity by 
virtue of its distance dependence. In fact, Ref. 2 shows a 
correlation function for the shadow fading that depends on 
distance. Thus, for both the fast and slow fading, the key to 
the averaging interval is vehicle velocity or its proportional 
surrogate, maximum Doppler frequency. Our approach to 
the adaptive averaging interval is, thus, through estimating 
the maximum Doppler frequency. 

The problem being discussed is to be distinguished from 
that addressed in [1] (also see Section 13.3.1 of [3]). In that 
reference, an averaging interval is being sought which will 
average out the effect of the Rayleigh fluctuations. By 
conOTist, we wish to exploit the Rayleigh fluctuations 
themselves to obtain information about the propagation 
environment and then use that information to determine the 
averaging interval. To further put our investigation into 
perspective, we are considering here an ideal Rayleigh fading 
environment. Further analysis will treat other realistic 
factors. One such factor is the distinction between signal 
strength and CA measurements (see [3], Section 9.1.2, in this 
regard). Furthermore, it should be recognized that other 
measures (e.g., bit error rate) are candidates for handover 
triggers (either in conjunction with, or instead of, signal 
strength). The present study should help in comparative 
evaluations of these measures by indicating how information 
inherent in the measures can be further exploited. 

There are a number of other considerations in handover 
decisions besides averaging intervals (see, e.g., [4]). It 
should also be observed that the three factors above need to 
be considered in power control measurements. In [5], the 
Rayleigh fluctuations are averaged out, while in [6], the 
DS/CDMA power control responds to the Rayleigh 
fluctuations fit needs to respond to relatively slow Rayleigh 
fading because the combination of interleaving and coding 
then becomes ineffective in smoothing out the fading 
fluctuations; see [6]). 
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The paper is organized as follows. In Section 2 the 
averaging interval tradeoff is discussed Section 3 shows a 
basis for choosing the averaging interval with information 
available in the measurements. Error estimates are given in 
Section 4. Adaptive use of this information is considered in 
Section 5. Comparison of our approach for estimating 
maximum Doppler frequency to some approaches in the 
literature is given in Section 6. Concluding remarks are given 
in Section 7. 

2. The Averaging Interval Tradeoff 

In [1], the standard deviation of time averaged Rayleigh 
faded signal strength measurements is derived. (Stationarity 
is assumed over the measurement intervals in question). This 
result is used in [1] to determine the minimum averaging 
interval for a given standard deviation of the measured power 
in dB . From inspection of Figure 2 of [1], the following 
approximation is found for the minimum averaging interval, 
T (k J d) seconds , for an averaged measurement standard 
deviation of k dB (i.e., an accuracy of k dB) and a 
maximum Doppler frequency of f d (in Hz): 

Remark (2-1). The maximum Doppler frequency, f dt 
equals velocity divided by wavelength, v/X. With T a time 
interval, /jT measures distance in number of wavelengths. 
While we shall be primarily stating results in terms of time, 
there is always the translation to the physically meaningful 
distance in number of wavelengths Thus, e.g., for an 
accuracy of 1 dB in (2-1), 20 wavelengths must be 
traveled. □ 

The required averaging interval T decreases with f d 
which is fortunately consistent with needing to make a 
quicker decision with faster vehicles. From (2-1), it also seen 
that, for a given standard deviation, the minimum averaging 
interval can be easily found and updated if /j can be easily 
evaluated and updated. 

In [1], the averaging interval tradeoff is affected by slow 
(lognormal) fading as well as fast (Rayleigh) fading. The 
specific example given in Figure 4 of [1] assumes a 0.2Hz 
maximum frequency for the spectrum of the slow fading. In 
[2] (Section 2), the bandwidth for the slow fading is given as 
proportional to vehicle velocity and, thus, proportional tof d . 
Hence, f d is the key parameter for the averaging interval 
taking account of both fast and slow fading. 

Remark (2-2). When averaging for both slow and fast 
fading, the slow fading predominates. □ 

Remark (2-3). The proportionality between f d and the 
key parameters for slow fading may have a direction- 
dependence, not yet investigated-^ 



There are various ways to estimate f d . Some possible 
approaches can be based on 

1. Mobile velocity (proportional iof d through the carrier 
frequency). 

2. Inferring f d from level crossings using results such as 
in f7], Chapter 6. 

3. Performing a spectrum or autocorrelation analysis and 
inferring f d from the spectrum or autocorrelation. 

4. Squared deviations of logarithmically compressed 
signal envelope (amplitude) measurements. 

We shall define and develop the last option in the next 
section. 

3. Using Squared Deviations of Signal Strength 
Measurements 

The rationale for using deviations of signal envelope 
amplitude measurements as in Point 4 above is that, if there 
is a fixed time interval between those measurements, it would 
be expected that the size of the deviations would increase as 
the fading becomes faster, i.e., as f d increases. We shall 
quantify that effect now. 

Remark (3-1). Note that we are using discrete-time 
measurements in contrast to the continuous-time 
measurements analyzed in [1]. For the sample mean 
examined in [1], the distinctions between discrete-time and 
continuous-time measurements are analyzed in a 
straightforward manner and usually slight and not otherwise 
consequential (see, e.g., p. 176 of [8]). By contrast, with our 
use of deviations here, the time difference between 
successive measurements will be seen to be key to the 
development of the method. Note that the averaging interval 
for our method consists of multiples of the time between 
successive (individual) measurements-^ 

Now, define 
where 

y, = 20 log l0 Xi (3-2) 

and xi is the i th signal envelope amplitude measurement 
(instantaneous, not averaged). Logarithmic compression 
facilitates estimation of f d since it is not necessary to make 
separate measurements of the mean (as will be seen in the 
following). For Rayleigh fading, the expectation of V is 
found to be approximately 

E ( V ) = 62.05 [ 1 -J 0 H2nf d x)\ (3-3) 

where t is the time between successive signal strength 
measurements (and, thus, the averaging interval is T~Nx). 
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For (3-3), we used the following approximation from [I]: 

cov(y i+k ji) = cov[y (t+kx)y(t)} = iS.S7^J^i2nf d kx) (3-4) 

It is thus seen that f d can be inferred from V. That is, 
/ d . the estimate of f d , is obtained from 

V= 62.05 [WoWrfT)] (3-5) 
Remark (3-2). (3-4) is actually an upper bound for the 
covariance function. It was found satisfactory for its use in 
[1]. For use in estimating f d , however, the estimate should 
be corrected based on numerical comparisons with the exact 
expression. The use of (3-4) (with the understanding that a 
correction will be made) simplifies the exposition here.D 



The inference of (3-5) is most transparent when 
2Kf d Z oc 1. 



(3-6) 



Then a series approximation for the Bessel function yields 

Even if (3-6) does not hold, it is easy to infer f d from V 
using (3-5) and we do not assume (3-6) in the analysis to 
follow. It should be recognized, however, that a condition 
Uke 27t/ rf x < 1 is needed for effectiveness of estimating /<, 
from (3-5) because of zeros of the Bessel function. This will 
be made evident from the error analysis of the next section. 
For/ d up to 100 Hz , a x of 1 ms yields 2nf d x < 0.628. 

4. Error Analysis 

An estimate for the the variance of V is (with 
*.=>»-£<>;)): 

var(V) = N- 1 var(y i+l - yi )2 

+ 2 N-S-pN-J)cov & M -n )2 (y.-^-y^) 2 ] 

=N-» variz^-zi)* 

+ 2 N~*"-pN-j )[cov (z l+1 2^ +>+1 2) 

-2C0V (z/ + i 2 ,Z£ +/ Z /+ ; +1 >+COV (Z/ + i 2 ,Z/ + ; 2 ) 

-2cov (z, z Ul ,Zi+ J+l *)+4cov (zi z Ux auj z UJ +1 ) 

-2COV {Zi Z U , ,Z £ +j 2)+COV (2; 1 2 ) 
-2C£>V (Zi l t z i+J Z i+j+1 }±COV (Zi 2 2)] 

+4 N-^'pN -j )(4 cov (y f y i+j )[cov (y, y i+j ) 

-cov (yi juj-d-cov (y,- o>/+y+i)J 
—+ [cov ft M+j^y+cov (y, y i+j + x y?) (4 -l) 

where i is arbitrary for a wide-sense stationary process and 
we have used (8-121) in [9]. This uses a joint Gaussian 
approximation for the z it as is commonly done in related 
correlation analyses (e.g., [10], p. 271). 

Remark (4-1). There is also a bias in the estimate of V 
due to nonlinearity which, however, is not as significant as 



the bias due to measurement noise. We shall return to this in 
Section 6.D 

An estimate for the relative error in f d , is given by (using 
(3-3) above and (5-56) of [1 1]): 

gC/jj - vW(V) 

fd 



(4-2) 



248.2*/,* \Jdi2%f d xyj x (2%f d x) \ 
Note, in particular, the unboundedness of this error 
estimate near zeros of the Bessel functions (Jq has a zero at 
2.4. J ! has zeroes at 0 and 3.8). Thus, one should choose x 
small enough to avoid those zeros. This confirms the 
comment about the size of t at the end of the last section. 
For the rest of this paper, we shall assume T = 1 ms . - — 

Numerical results for this relative error are given in 

Figure 1. The relative error, S^kl, is plotted versus the 

J d 

averaging interval, T (in seconds), with f d (in Hz) as a 
parameter, t is fixed at 1 w (for the reason given in the last 
paragraph). The use of (3-7) facilitates interpretation of the 
results. From (3-7), we obtain 



o(fd) 
E(fd) 



<j(jV) 
£(VV) 



(4-3) 



(4-4) 



and, then -relating- the* coefficient of variation of a random 
variable to that of its square root, 

g(A<) - 1 o(V) 
E(f d ) ~ T E(V) 
Remark (4-2). The errors in Figure 1 are to be taken as 
qualitative indicators. First of all, several approximations 
were used to derive them (simulations indicate that they are 
underestimates). Secondly, measurement noise needs to be 
considered (see Remark (4-l)).D 

5. Adaptive Measurements 

The results of Section 3 indicate that the averaging 
interval can be done adaptively. In particular, we would like 
to adaptively choose the number of samples N (with x fixed) 
for the sample mean S 



(5-1) 



Remark (5-1). The continuous-time version of Eq. (5-1) 
is examined in [1]; see Remark (3- 1).D 

One convenient way to facilitate that is to replace (5-1) 
by exponential smoothing (e.g.. [12], Section 4.3) and 
continually adjust the weights: 

Si+i = a(V.>y i+ i + [l-aWtfi (5-2) 
= P(Vi Xfc+i-yo 2 + [1-VtViWi (5-3) 

Here, a and 0 are both exponential weighting factors, each 
continually updated by the most recent V estimate, which is 
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related to/j. The functional forms of these weighting factors 
are to be determined for differing situations. A general basis 
for the choice is as follows. 

1 . Given Vj , estimate ft . 

2. From/j, determine an appropriate N for (5-1). This 
involves the tradeoff between estimation error and lag 
in the h an doff algorithm response, as discussed. This 
tradeoff is not resolved here, requiring more work 
which also accounts for the combined roles of 
averaging interval and hysteresis (see [13], [14]). 

3. Relate a(v;) to iV(v;). This is done by noting that, for 
constant a, (5-2) is equivalent to 



(5-4) 



so that the past sample y;_y is attenuated by (l-a)o/. 
By setting (l-oOa'' equal to say, lie, one can 
determine an an a corresponding to an N . 

4. Use aiVi ) to update (5-2). 

5. Determine p(V; ) in a similar manner to update (5-3). 

6. Relationship to Some Other Work on Doppler Estimation 

In [15], techniques are given for estimating the Doppler 
spread parameter D, defined as 



4 



J(v-v) 2 />(v)dv 



}P(v)dv 

where P (v) is the power density spectrum (with respect to 
time) of T(f ,t), the received narrowband process in response 
to a complex sinusoid. This is a general definition of 
Doppler spread while we are using the specific form of 
maximum Doppler frequency encountered in mobile radio. 
Estimates of D in [15] are given in terms of derivatives. In 
actually calculating derivatives, finite differences would 
usually be used. This is explicit in [16] (specifically, Eq. 
(24)), where further considerations are given for estimating 
O. By contrast, we do not use differences as approximations 
for derivatives but directly exploit the property of the 
differences. 

The accuracy of one of the approaches of [15] is 
considered in [17]. They consider three sources of error: 

1. Noise in measurements. 

2. Sampling error. 

3. Finite difference approximation to derivative. 
The relationship to our work is: 

I. This results in a significant bias which must be 
incorporated into the present study. 



2. Independent samples are considered in [17] while we 
analyzed correlated samples. 

3. We do not try to approximate a derivative with finite 
differences, as mentioned above. 

7. Concluding Remarks 

We have explored the derivation of information from the 
signal strength variations that is useful for adaptive handover 
algorithms. Further work is needed to test the idealizations 
assumed and to use the approach discussed here (and/or 
others such as mentioned at the end of Section 2) in a real- 
time handover algorithm. The idealizations include a purely 
Rayleigh environment with noiseless measurements. A 
number of sensitivity studies are needed. Furthermore, the 
interactions between averaging interval and hysteresis must 
be considered (see [13], [14]). It will also be of interest to 
see how a modification of the algorithm works in an 
environment dominated by non-Rayleigh effects such as 
turning corners. 

As mentioned in Section 2, can also be obtained from 
vehicle velocity if that information is available to the handoff 
algorithm processing system. On the other hand, an estimate 
of fd might be used to estimate velocity. 
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Figure 1 - Relative Errors 
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